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We describe the effective heavy meson lagrangian for S- and P-wave heavy-light mesons in terms 
of a model based on meson-quark interactions, where mesonic transition amplitudes are represented 
by diagrams with heavy mesons attached to loops containing heavy and light constituent quarks. 



I. INTRODUCTION 
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The model is relativistic and incorporates the heavy quark symmetries. The universal form factors 
of the heavy meson transition amplitudes are calculated together with their slopes and compared 
to existing data and limits. As further applications of the model, strong and radiative decays of D* 
and B* are considered. The agreement with data is surprisingly good and shows that the model 
offers a viable alternative to effective meson lagrangians which require a larger number of input 
parameters. 
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An effective theory for heavy mesons, implementing the heavy quark symmetries, has been very successful at the 
phenomenological level (for reviews see jjj and references therein). Predictions are easily obtained once the unknown 
ij effective couplings are fixed from experimental data. Moreover such an effective approach can be combined with chiral 
symmetry for light mesons, thus giving a simple framework for implementing the known approximate symmetries of 
\ quantum chromo dynamics (QCD) [||. The disadvantage of such an approach is the number of free parameters, which 
grows very rapidly if one tries to improve the calculations beyond leading order. 

In order to go beyond the symmetry approach, one should be able to derive the free couplings at the meson 
lagrangian level from a more fundamental theory, for example directly from QCD. This is clearly a difficult task. In 
the long run the definitive answer will come from systematic first principle calculations, such as in lattice QCD, but 
at present heavy meson physical quantities such as spectra and form factors are still subject to extrapolations, even if 
recent improvements are impressive g]. Moreover an alternative and intuitive way for interpolating between QCD and 
an effective theory would be interesting in itself, allowing us to understand better the underlying physics. Obviously 
there is a price to pay for any simplification that may allow to calculate the parameters of an effective heavy meson 
theory without solving the nonperturbative QCD problem. Our point of view here will be to consider a quark-meson 
lagrangian where transition amplitudes are represented by diagrams with heavy mesons attached to loops containing 
heavy and light constituent quarks. It should be kept in mind that what we study here is only a model and not full 
QCD. However one can hope to describe the essential part of the QCD behavior in some energy range and extract 
useful information from it. 

Since the model used in the present paper is based on an effective constituent quark-meson lagrangian containing 
both light and heavy degrees of freedom it is constrained by the known symmetries of QCD in the limit ttiq — > oo. 
A similar model, for light quarks only, was pioneered, in the context of chirally symmetric effective theories, several 
years ago in ||. A related approach is the one based on the extended Nambu Jona-Lasinio (NJL) models ||, 
whose generalization to include both heavy and light quarks has been studied in B. In this approach path integral 
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bosonization is used, which replaces the effective 4-quark interactions by Yukawa-type couplings of heavy and light 
quarks with heavy meson fields. One may think of going even further and try to bosonize directly QCD, but a number 
of approximations are needed and one usually ends up with a non-local lagrangian. 

In the following we are interested in the description of the heavy meson interactions for the doublets with spin/parity 
J p = (l + ,0 + ), (1 _ ,0~), (2 + , 1 + ), i.e. S- and P-wave heavy-light mesons. These states are the object of experimental 
search Q. Our aim is therefore to extend the approach of [Q and |^ to the description of the three spin-parity 
doublets mentioned above. We will consider a lagrangian at the meson-quark level, fixing the free parameters from 
data. This will allow us to deduce from a small number of parameters the heavy meson couplings and form factors, 
with a considerable reduction in the number of free parameters with respect to the effective lagrangian written in 
terms of meson fields only. 

The plan of the paper is as follows. In section 2 we give an outline of the model. In Section 3 we compute some 
fundamental parameters of the model, such as mass splittings and renormalization constants. The model is used in 
Section 4 to obtain predictions for semileptonic B decays into negative parity and positive parity charmed resonances; 
in particular we compute the relevant universal form factors that describe such decays in the uiq — > oo limit and 
discuss the Bjorken sum rule ||. In section 5 we compute strong and radiative D* and B* decays. Finally in Section 
6 we draw our conclusions. 



II. THE MODEL 



To begin with we briefly review the description of the heavy quark field in the heavy quark effective theory (HQET) 
JD. The heavy quark is defined as follows: we indicate by i> M the velocity of the hadron containing the heavy quark 
Q. This is almost on shell and its momentum pq can be written, introducing a residual momentum fc, as 

PQ = m Q v + k. (1) 

One extracts the dominant part rriQV of the heavy quark momentum by defining the new field Q v 

Q v (x) — e-xjp(irriQVx)Q(x) — h v (x) + H v (x). (2) 

The field h v is the large component field, satisfying the constraint fh v ~ h v . If the quark Q is exactly on shell, it 
is the only term present in the previous equation. H v , the small component field, is of the order 1/toq and satisfies 
^H v = —H v . The effective lagrangian can be derived by integrating out the heavy fields in the QCD generating 
functional pj, however here we only show the effect of the niQ — > oo limit on the Feynman rules relevant to our 
calculation . The heavy quark propagator in the large itiq limit is: 

' ~ 1 (3) 



f — m,Q v ■ k 2 

The (1 + ^)/2 projection operator can always be moved close to a spinor h v satisfying fh v = h v since f commutes 
with the heavy quark-gluon vertex. 

In our model the heavy quark propagator appears in the loops always together with light quark propagators as the 
heavy mesons we consider are made up of constituent heavy quark and light antiquark. The light-quark momentum 
k may be taken equal to the integrated loop momentum and the heavy-quark propagator in the heavy quark limit 
becomes: 

(4) 



v ■ k + A 

where A is the difference between the quark mass and the common mass of the heavy meson doublet. 



A. Heavy meson fields 



In order to implement the heavy quark symmetries in the spectrum of physical states the wave function of a 
heavy meson has to be independent of the heavy quark flavor and spin. It can be characterized by the total angular 
momentum sg of the light degrees of freedom. To each value of S£ corresponds a degenerate doublet of states with 
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angular momentum J = Sf, ± 1/2. The mesons P and P* form the spin-symmetry doublet corresponding to sg = 1/2 
(for charm for instance, they correspond to D and D*). 

The negative parity spin doublet (P,P*) can be represented by a 4 x 4 Dirac matrix H, with one spinor index for 
the heavy quark and the other for the light degrees of freedom. These wave functions transform under a Lorentz 
transformation A as 



H -> D(A) P D(A) _1 



(5) 



where -D(A) is the usual 4x4 representation of the Lorentz group. Under a heavy quark spin transformation S, 
H — > S*iJ , where S satisfies [jtf, S] = to preserve the constraint /P = iJ. 
An explicit matrix representation is: 



H = 



(1 



^7s 



Here w is the heavy meson velocity, v^P* = and Mh = Mp = Mp- . Moreover -^H = — Hj/f = H, Htf 
and P*^ and P are annihilation operators normalized as follows: 

<0|P|Qg((T)) = ^fm I 



|P*1Q?(1-)) 



'A/ 



(G) 
(7) 

-0 = H 



(8) 
(9) 



The formalism for higher spin states is given in [ p"T| . We will consider only the extension to P- waves of the system 
Qq. The heavy quark effective theory predicts two distinct multiplets, one containing a + and a 1 + degenerate state, 
and the other one a 1 + and a 2 + state. In matrix notations, analogous to the ones used for the negative parity states, 
they are described by 



S 



1+^ 



[Pfrf'T* - Po] 



(10) 



and 



r 2 iv 



These two multiplets have se = 1/2 and sg — 3/2 respectively, where sg is conserved together with the spin sq in the 
infinite quark mass limit because J = sg + sq. 

B. Interaction terms 



The effective quark-meson lagrangian we introduce contains two terms: 

C = Cgg + Chi ■ 



(12) 



The first term involves only the light degrees of freedom, i.e. the light quark fields \ an d the pseudo-scalar SU (3) 
octet of mesons n. At the lowest order one has: 



Igg = x(iD^ + 5A^ t 7 Al 75)x - m XX 
f 2 



(13) 



Here D M = 9 M - iV„, £ = exp(wr//^), E = £ 2 , U = 130 MeV and 



V" = ^(^^C + ^t) 



(14) 
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The term with qa is the coupling of pions to light quarks; it will not be used in the sequel. It is a free parameter, but 
in NJL model $a — 1- Differently from Q the gluon field is absent in the present description. Apart from the mass 
term Cu is chiral invariant. 

We can introduce a quark-meson effective lagrangian involving heavy and light quarks and heavy mesons. At lowest 
order we have: 



Cfu = Qviv ■ dQ v - \ x{H + S + iT ti — )Q v + h.c 

+ 2^ Tr ™ + 2^ Tr[55] + 2^ Tr[ ^ TAl] (15) 

where the meson fields H, S, T have been defined before, Q v is the effective heavy quark field, G3, G' 3 , G4 are 
coupling constants (to allow a comparison with previous work we follow, as far as possible, the notations of ||) and 
A x (= 1 GeV) has been introduced for dimensional reasons. Lagrangian (15) has heavy spin and flavor symmetry. 
This lagrangian comprises three terms containing respectively H , S and T. From a theoretical point of view, the G3 
and G 3 couplings could have been determined by Ah and As, the experimental values of the mass differences between 
the mesons, belonging to the H and S multiplets respectively, and their heavy quark constituent masses. Since these 
data are not all known, we have adopted the dynamical information coming from the NJL-model as proposed in j(| 
and so we use the following modified version of ( |l5|) in which we give the fields H and 5* the same coupling to the 
quarks: 

/ _ _ D* 1 

Cm = Qviv ■ dQ v - \ r(H + S + iT ll —)Q v + h 

+ 2^ Tr t( 5 + S)(H- S)] + ±-Tv[T^] (16) 

In doing this we assume that this effective quark meson lagrangian can be justified as a remnant of a four quark 
interaction of the NJL type by partial bosonization. H and S are degenerate in mass in the light-sector chirally 
symmetric phase. In the broken phase the mass splitting is calculable and will be computed in the next section 
in terms of the order parameter m, the constituent light quark mass. The part containing the T field can not be 
extracted from a bosonized NJL contact interaction and requires a new coupling constant G4. 

An essential ingredient of the nonperturbative behavior of QCD is the suppression of large momentum flows through 
light-quark lines in the loops. The way this suppression is introduced in the model is crucial and is part of the definition 
of the model itself. One can for example include factors in the vertices which damp the loop integrals when the light- 
quark momentum is larger than some scale of the order of 1 GeV |l^] . Our model describes the interactions in terms 
of effective vertices between a light quark, a heavy quark and a heavy meson; we describe the heavy quarks and heavy 
mesons consistently with HQET, and thus the heavy quark propagator in the loop contains the residual momentum 
k which arises from the interaction with the light degrees of freedom. It is therefore natural to assume an ultraviolet 
cut-off on the loop momentum of the order of A ~ 1 GeV, even if the heavy quark mass is larger than the cut-off. 

Concerning the infrared behavior, the model is not confining and thus its range of validity can not be extended 
below energies of the order of Aqcd- In practice one introduces an infrared cut-off /i, in order to drop the unknown 
confinement part of the quark interaction. 

We choose a proper time regularization (a different choice is followed in |l3|). After continuation to the Euclidean 
space it reads, for the light quark propagator: 

d 4 fc g >2 1 = - / d 4 k E I' 1 * dse~ s ^+ m ^ (17) 
k 2 E + m 2 J 

where fi and A are infrared and ultraviolet cut-offs. This choice of regularization is consistent with the one in the 
second paper in B. In (17) m is the constituent light quark mass, whose non vanishing value implies the mass splitting 
between the ff-type and S'-type multiplets. Its value, as well as the values of the cutoffs /i and A are in principle 
adjustable parameters. On physical grounds, we expect that A ~ 1 GeV, m ~ fi ~ 10 2 MeV. Given our choice of the 
cut-off prescription as in we assume for A the same value used there i.e. A = 1.25 GeV; we observe, however, that 
the numerical outcome of the subsequent calculations are not strongly dependent on the value of A. The constituent 
mass m in the NJL models represents the order parameter discriminating between the phases of broken and unbroken 
chiral symmetry and can be fixed by solving a gap equation, which gives m as a function of the scale mass fi for given 
values of the other parameters. In the second paper of ref. the values m = 300 MeV and fi = 300 MeV are used 
and we shall assume the same values. As shown there, for smaller values of fi, m is constant (=300 MeV) while for 
much larger values of /it, it decreases and in particular it vanishes for fj, = 550 MeV. 
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III. RENORMALIZATION CONSTANTS AND MASSES 



For the computation of the constants G3, G4 appearing in the quark-meson effective lagrangian of Eq.(15) we 
consider the meson self energy diagram depicted in Fig.l. 



Fig. 1 - Self-energy diagram for the heavy meson field. 

For H and S states the evaluation of Fig.l gives the results (k is the residual momentum, see eq. (Q)): 

d 4 £ Tt[{/-^-m)HH] 



Ti[-HU H (v ■ k) H) = iN c 



and 



Tr[+S n s (v ■ k) S] = iN c 
respectively, and for T states we get 



(2tt) 4 {{£ - k) 2 -m 2 + ie)(v£ + ie) 

d A £ Tr[(/ - A? + m)SS] 

(2tt) 4 ((£ - k) 2 ~m 2 + ie)(v£ + ie) 



Tr[+T„ U T {v ■ k) T p ] = 



iN c f dH Tr[(£ M - fc M ) 



J + m)(t v - k v )T»T v ] 



K J (27T) 4 



((£ - k) 2 - m 2 + ie)(v£ + ie) 



(18) 



(19) 



(20) 



We notice that the only difference between Hh and II5 is in the sign in front of m in the numerators of the 
integrands. We now expand (p~8|) , (p~S|) andfl2(i|) around A#, Ag and At respectively: 



IL(v ■ k) ~ n(A) + n'(A)(« • k - A) . 
In this way we obtain the kinetic part of the effective lagrangian for heavy mesons in the usual form Q 

C ef / = TtH(iv -d-A H )H + TrS(iv -d-A s )S 



(21) 



TiT^iv ■ d- A T )T^ 



provided we satisfy the conditions: 



1 
1 

2G7 



n H {A H ) = n s (A s ) 
n T (A T ) 



and renormalize the fields as 



H 

S 
f 



H 

s 



Z s 
T 



' Zt 



The renormalization constants are given by: 



(22) 

(23) 
(24) 

(25) 
(26) 
(27) 

(28) 



Here, prime denotes differentiation and j = H, S and T. In the previous equations, Ah, As, At are the mass 
differences between the heavy mesons H, S, T and the heavy quark. The expressions of n#, II5, and Ht are: 
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U H (A H ) =h + (A H + m)I 3 (A H ) 
U s (As) = h + (A s -m)h(A s ) 



n T (A T ) = — 



■[I (A T ) + A T h + (A T 



'■)h(A 7 



while the field renormalization constants are 



Z 



H 



(A H + m) 



dh(A H ) 
OA 



H 



h(A H ) 



Zs X = (As-m)^M + I 3 (A s ) 



dA s 



(29) 
(30) 

(31) 



(32) 



(33) 



ZZ 



3AI 



+ (m - 



A T )[ 



- m 2 )[(m 

dI (A T ) 
dA T 



A T ) 



h 



dh (At) 



dA T 
2A T h(A T )] 



h (A T )] 



A T h 



(34) 



The integrals Iq,Ii,I[ and J3 are defined in the Appendix. We employ the first equation in ( |23| ) to obtain G3 from a 
given value of Ah (we use Ah in the range 0.3 — 0.5 GeV), while the second equation in ( p3| ) can be used to determine 
A s . 

For At we can use experimental information. Let us call Mh and Mt the spin-averaged ii-multiplet and T-multiplet 
masses. We can write Mh = tuq + Ah + A' H /mQ and Mt = tuq + At + A t /itiq. For charm, the positive parity 
experimental masses and widths are as follows: for we have mr>*(2460)° — 2458.9 ± 2.0 MeV, Tp*^^) = 23 ± 5 



MeV and to D j(2460)J 
belonging to the T § 



= 2459 ± 4 MeV, r D »( 2 4 6 o)± = 25+* MeV. These particle are identified with the state 2+ 
multiplet. As for the Z)* (2420) state, experimentally we have m£>*(242o)° = 2422.2 ± 1.8 MeV, 
T D * (2420)° = 18.9+3 5 MeV; this particle can be identified with the state 1 + belonging to the T multiplet, neglecting 
a possible small mixing between the two 1 + states belonging to the S and T multiplets We observe that the 
narrowness of the states D% and D\ is due to the fact that their strong decays occur via D-waves. From this analysis 
we get At — Ah + (A' T — A' H )/m c ~ 470 MeV. For beauty, the experimental data on positive parity resonances show 
a bunch of states with a mass M B *> = 5698 ± 12 MeV and width r = 128 ± 18 MeV @. If we identify this mass 
with the mass of the narrow T states, we obtain At — Ah + (A' T — A' H )/mb ~ 380 MeV. For reasonable values of 
the heavy quark masses, we get 

At — A H ^ 335 MeV , (35) 
which is the value we adopt. The preceding analysis produces the figures in Table I. 



A H 


A s 


At 


0.3 


0.545 


0.635 


0.4 


0.590 


0.735 


0.5 


0.641 


0.835 



TABLE I. A values (in GeV) 



In Table II wc report the computed values of Gj and Zj for three values of Ah- 

Let us finally observe that we predict, according to the value of Ah, a value for the S— multiplet mass m = 2165±50 
MeV; these states, called in the literature Dq,DI' have not been observed yet since they are expected to be rather 
broad (for instance in |L5| one predicts r(fJ> -> D+tt") ~ 180 MeV and T(Dl' -> D*+n-) ~ 165 MeV). Theoretical 
predictions in the literature are somehow larger (m ~ 2350 MeV). 



A H 


I/G3 


Z H 


Z s 


Zj 1 


I/G4 


0.3 


0.16 


4.17 


1.84 


2.95 


0.15 


0.4 


0.22 


2.36 


1.14 


1.07 


0.26 


0.5 


0.345 


1.142 


0.63 


0.27 


0.66 



TABLE II. Renormalization constants and couplings. Ah in GeV; G3, Ga in GeV 2 , Zj in GeV 
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IV. SEMILEPTONIC DECAYS AND FORM FACTORS 



To begin with, we compute the leptonic decay constants F and F + that are defined as follows: 



(0\q^j 5 Q\H{0-,v)) = i^Mhv^F 
{Q\q^Q\S{Q + ,v)) = iy^v^F + 



(36) 
(37) 



It is easy to compute these constants by a loop calculation similar to that considered in the previous section for the 
self-energy ; one finds 



F : 

F+ 



G 3 
G 3 



(38) 
(39) 



The numerical results for different values of the parameters can be found in Table III. 





F 


F+ 


0.3 


0.33 


0.22 


0.4 


0.34 


0.24 


0.5 


0.37 


0.27 



TABLE III. F and F + for various values of Ah- A h in GeV, leptonic constants in GeV 3/2 
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We observe that, neglecting logarithmic corrections, F and F + are related, in the infinite heavy quark mass limit, 
to the leptonic decay constant f B and / + defined by 

(0\q^ 5 b\B(p)) = ip^f B (40) 

(0\q^b\B o (p)) = ipff+ (41) 

by the relations f B = F/y/M^ and /+ = F + /y/M Bo . For example, for A H = 400 MeV, we obtain from Tabic III: 

f B ~ 150 MeV (42) 

/+ ~ 100 MeV . (43) 

The numerical values that can be found in the literature agree with our results. For example the QCD sum 
rule analysis of |l(| gives F — 0.30 ± 0.05 GeV 3 / 2 (without a s corrections) and higher values (around 0.4 — 0.5 
GeV 3 / 2 ) including radiative corrections. As for lattice QCD, in || the value summarizing the present status of lattice 
calculations for f B is 170 ± 35 MeV. As for F+, a QCD sum rule analysis Q gives F + = 0.46 ± 0.06 GeV 3 / 2 , which 
is significantly higher than the results reported in Table III. 

Let us now consider semileptonic decays. The first analysis to be performed is the study of the Isgur-Wise function 
£ which is defined by: 

(D(v')\c^(l - -Ys)b\B(v)) = y/M B M D C cb + <) (44) 

(D*(v',e)\c^(l - ~f 5 )b\B(v)) = y/M B M D .Ccb ^)[i^ va0 e* u v' a v p 

-(l+ w )e*+(e*.tO«;j (45) 

where u) = v ■ v' and C c b is a coefficient containing logarithmic corrections depending on a s ; within our approximation 
it can be put equal to 1: C c b = 1. We also note that, in the leading order we are considering here £(1) = 1. To 
compute £ one has to evaluate the diagram of Fig. 2. 
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The £ function is plotted in Fig. 3. 

The integrals I3, I5 can be found in the Appendix. It is obvious from (|32|) , ([f6|) that £(u>) is correctly normalized, 
i.e. £(1) = 1. For the slope of the Isgur-Wise function, defined by 



P 2 IW = -|(D, (47) 

one gets the values reported in Table IV. The plot of the Isgur-Wise function is given in fig. 3 for three values 
of A#. It is useful to compare our results with other approaches. QCD sum rules calculations are in the range 
p 2 w = 0.54 — 1.0, therefore they agree with our results pjjfl . The results from different quark models ]l9| , po[ , 
p2[ l have been recently reviewed by the authors of [^3| ; they have studied the results of different models with the 
common approach of boosting the wave- functions by the Bakamjian-Thomas method [ p4[ ; for p 2 w they obtain results 
in the range 0.97— 1.28, to be compared to the results in table 4. Lattice QCD gives significantly smaller results; for 
example in |2^] the value p 2 w (l) = 0.64 is obtained. 

Let us now turn to the form factors describing the semileptonic decays of a meson belonging to the fundamental 
negative parity multiplet H into the positive parity mesons in the S and T multiplets. Examples of these decays are 

B -> D**£u (48) 

where D** can be either a S state (i.e. a + or 1 + charmed meson having sg — 1/2) or a T state ( i.e. a 2 + or 1 + 
charmed meson having si — 3/2). 

The decays in ([l8|) are described by two form factors 1^/2, t 3 / 2 (2(| according to 

< D* 2 (v', e)|c 7M (l - j 5 )b\B(v) >= yjZM B M D . t 3/2 (uj) x 

[ie^e^v'^ - [(w + l)e* m v a - (49) 



M B M D * 

<D* 1 (v', e)\^(l- 75 )b\B(v) >= ^ ^ r 3/2 H x 

{(oj 2 - l)e* + (e* • w)[3 v M - (oj - 2)t£] - i(w + lJe^e'W} (50) 



(A,(u')|c7 M (l ~ 75)fc|B(«)) = y/M B M Do 2 r 1/a («)(t^ - w M ) 
(£>IV, e)|c 7p (l - 75)6|S(w)) = yj ' M B M DV r 1/2 (w){2 i e^e* W 
+2[(1- W )e' + (e* •«)<]} . 



(51) 



(52) 



In all these equations we neglect logarithmic corrections. The form factors Ti/ 2 (oj), t 3 / 2 (uj) can be computed by a 
loop calculation similar to the one used to obtain £(oj). The result is 



n/aM = [h(A s ) - h(A H ) + (A H - A s + m(l - u)) h(A H ,A s , u>)] 

2(1 — OJ) 



and 



73/2(^) 



V^h^t r fI 3 (A H )~I 3 (A T )-{A H ~A T )I 5 (A H ,AT,oj) 
m\ 



y/3 



2 (1-oj) 



h(A H ) + h(A T ) + {A H + A T ) I 5 (A h ,A t ,oj) 



1 



2 (1+w) 

3S(A h ,A t ,oj)-(1-2uj) S(A t ,A Hi lu) 



2 (-1 - uj + uj 2 +u 3 ) 
+ (1 - uj 2 ) T(A H , A t , uj) - 2 (1 - 2u) f/(A H , A T , w 



where the integrals S 1 , T, {/ are defined in the appendix. 

Given the small phase space which is available for these decays (lu„ 
D*', Do), we can approximate 



(53) 



(54) 



1.33 for D*, D* 2 and u max ~ 1.215 for 
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r^-r^x [l-pf( W -l)] 



(55) 



Numerically we find the results reported in Table IV. 







Piw 


n/a(i) 


Pl/2 


T 3/2(l) 


P3/2 


0.3 


1 


0.72 


0.08 


0.8 


0.48 


1.4 


0.4 


1 


0.87 


0.09 


1.1 


0.56 


2.3 


0.5 


1 


1.14 


0.09 


2.7 


0.67 


3.0 



TABLE IV. Form factors and slopes. A H in GeV. 



An important test of our approach is represented by the Bjorken sum rule, which states 



2 

Piw 



(56) 



Numerically we find that the first excited resonances, i.e. the S and T states (k = 0) practically saturate the sum 
rule for all the three values of Ah- 

In Table V we compare our results (for Ah = 0.4 GeV) with other approaches. For r 3 / 2 we find a broad agreement 
with some of the constituent quark model results, whereas for ti/% we only agree with 



n/a(i) 



0.09 

0.41 

0.25 

0.31 

0.41 

0.059 

0.225 



~pj^7 



73/2(1) 



1.1 
1.0 

0.4 

2.8 

1.4 

0.73 

0.83 



0.56 

0Al(input) 

0.28 

0.31 

0.66 

0.515 

0.54 



P3/2 



Ref. 



2.3 

1.5 

0.9 

2.8 

1.9 

1.45 

1.50 



This work 
[27j 

\n\ 

[22] 
[28] 

[23], |l| 
[23], [19| 



TABLE V. Parameters of the form factors T1/2, t 3 / 2 . The results in this table are for Ah = 0.4 GeV. 

Finally in Table VI we present our results for the branching ratios of B semileptonic decays to S— and P— wave 
charmed mesons for three values of Ah computed with V c b = 0.038 |29|] and tb = 1.62 psec. We see that data favor 
a value of A H ^ 400 - 500 MeV. 



Decay mode 


Ah = 0.3 


A H = 0.4 


A H = 0.5 


B -> Dlv 


3.0 


2.7 


2.2 


B -> D*iv 


7.6 


6.9 


5.9 


B -> D £u 


0.03 


0.005 


0.003 


B -> D\'lv 


0.03 


0.008 


0.0045 


B -> D*£i/ 


0.27 


0.18 


0.13 


B -> D\lv 


0.43 


0.34 


0.30 



Exp. 



1.9 ±0.5 flj 
4.68 ±0.25 



0.74 ±0.16 
< 0.85 



TABLE VI. Branching ratios (%) for semileptonic B decays. Theoretical predictions for three values of Ah and experimental 
results (for B° decays). Units of Ah in GeV. 
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V. STRONG AND RADIATIVE DECAYS OF HEAVY MESONS 



In this section we consider the strong decays 



#-> Hit (57) 
Hi: (58) 



as well as the radiative decay 

P* -> P7 (59) 
where P* and P are the 1~ and 0~ members of the multiplet H. 

A. Strong decays 

The calculation of the strong coupling constants describing the decays H — > Htt (i.e. D* — > -D71") and 5 — * i?7r is 
strongly simplified by adopting the soft pion limit. In the case of the decay D* — > Z)7r this procedure introduces only 
a small error since the phase space is actually very small (for the S — > Hn decay the situation is different, see below). 

Let us define gD*D-K by the equation: 

< TT + (q)D°(p)\D* + (p', e) >= i<to* D7r e% . (60) 
The constant gwDir is related to the strong coupling constant of the effective meson field theory g appearing in 0] 

C = igTriHH^A^ + [ih TriHS^ysA^ + h.c] (61) 

by the relation 

2mD (RO\ 

gD'DTr = —^9 ( 62 ) 

valid in the mg — > 00 limit. 

To compute g in the soft-pion-limit (q^ — » 0) we consider the matrix element of d^A^ and derive a Goldberger- 
Treiman relation, following the approach of this approach differs from the method employed in the first paper in 
H, that assumes the so-called low-momentum-expansion approximation and a mixing between the Nambu-Goldstone 
bosons and higher mass resonances. Our method amounts to a loop calculation involving a current and two H states 
(see Fig. 4). 



Fig. 4 - Pion vertex on the light quark line for the heavy-heavy-pion interaction 
The result for g is as follows: 

.9 = Zh 



i/ 3 (A ff ) - 2(m + ^A ff )(7 2 + A H h(A H )) - ^m 2 h(A H ) 



where I2, I3 have been defined already and I4 is in the Appendix. Numerically we get 

g = 0.456 ± 0.040 



(63) 



(64) 
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where the central value corresponds to Ah = 0.4 GeV and the lower (resp. higher) value corresponds to Ah = 0.3 
GeV (resp. Ah — 0.5 GeV). These values agree with QCD sum rules calculations, that give g = 0.44 ± 0.16 (for a 
review see 0|; see also [j32|), with the result of relativistic constituent quark model: g ~ 0.40 31 1, g — 0.34 Q. From 
the computed values of g we can derive the hadronic width using 



*+ 



D°n+) 



9 



6VT/2 



\Pn\ 



(65) 



The numerical results will be discussed at the end of this Section. Let us now comment on the strong decay S — > Htt. 
A complete calculation of the corresponding decay constant h appearing in (61) is much more involved and one can 
preliminary try to compute it in the soft-pion limit. This approximation implies that we must assume As = Ah] 
putting A = (Ah + As)/2, we obtain 



h = Z{I 3 (A) +2A/ 2 + 2(A 2 



were Z is given by 



)4(A)} 



(66) 



Z = 



'3(A) + A— -Irn— 



We obtain, for A in the range 0.43 — 0.57 GeV the result 

h= -0.85 ±0.02 



(67) 



(68) 



which is somehow higher, but still compatible, within the theoretical uncertainties, with a result obtained by QCD 
sum rules: h = -0.52 ± 0.17 Jill. 



B. Radiative decays 

Let us now consider the radiative decays D* — > and B* — > Bj. To be definite we consider only the former decay; 
moreover we use the SU(3) flavor symmetry for the light quarks. The matrix element for this radiative transition is 
the following: 

M (D* -> Dj) = i en v „ p' aPl3 (69) 

where e M is the photon polarization and the coupling [i comprises two terms: 

H = n' + ti h , (70) 

corresponding to the decomposition: 

M(D* -^D 1 ) = e e*" < D(p')\je m \D*(p, n ) > 

= ee*» <D(p')\4 + J^D*{p,rj)> . (71) 

Here and are the light and the heavy quark parts of the electro-magnetic current: 

2 1 - 1 



and 



Correspondingly, eq. (70) becomes 



4 = - h ltL b = J2 cqQIhQ- (73) 

'* Q=c,b 



M = / + M' l -f + (74) 
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where A g and Aq are mass parameters to be determined. The photon insertion on the heavy quark line generates the 
fx coupling, while the /i e arises when the photon is inserted on the light quark line of the loop. 

In the itlq — > oo limit, the matrix element of can be expressed in terms of the Isgur-Wise form factor £(u>) as 
follows: 

< Dtf)\J*\D*fan) >= e c < D(p')\cY l c\D*(p, V ) > 
2 



- i-^M D M D ^(to)e^ a ^v' a v^ (75) 

where p' = Mjjv' , p — Me>*v and lo ~ 1 because: 

= q 2 = m% + M 2 D , - 2M D M D »v ■ v'. (76) 

Taking into account the normalization £(1) = 1 one gets, for the charm and beauty mesons respectively (with 
q = d), 

-i: - - 55 (77) 

with 



A c = y/M D M D * A b - y/M B M B , . (78) 

In the leading order in \/m c and I /nib one finds: 

A c = m c A b = m b (79) 

which fixes . As for the coupling in our approach it stems from the diagram where the photon line is inserted 
on the light quark propagator (see Fig. 5). 




Fig. 5 - Electro-magnetic current insertion on the light quark line. 
The result is: 

= A" 1 = 2 x Z H [I 2 + (m + A H )I 4 (A H )} (80) 

Numerically we get: 

(3 = 1.6 ±0.1 GoV" 1 (81) 

where the central value corresponds to Ah = 0.4 GeV and the lower (resp. higher) value corresponds to Ah = 0.5 
GeV (resp. A# = 0.3 GeV). These values agree with the result of Heavy Quark Effective Theory and the Vector 
Meson Dominance hypothesis j34| . 
One can use the formula: 

T(D* ^ D 7 ) = ^^\rf\k\* (82) 
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(k = photon momentum), to compute the D* radiative widths. Using this equation together with ( |65| ) we obtain the 
results of Table VII. 



Decay mode 


A H = 0.4 GeV 


A H = 0.5 GcV 


Exp. 


D* u -> D°TT° 


65.5 


70.1 


61.9 ±2.9 


D*° -> L>°7 


34.5 


29.9 


38.1 ±2.9 


D *+ £,0^ + 


71.6 


71.7 


68.3 ± 1.4 


D *+ 


28.0 


28.1 


30.6 ±2.5 


D *+ _> D + 1 


0.4 


0.24 


1 1+ 2 - 1 

L - l -0.7 



TABLE VII. Theoretical and experimental D* branching ratios (%). Theoretical values are computed with Ah = 0.4, 0.5 
GeV. 

Taking into account the approximations involved in the present calculation, we find comparison between theoretical 
prediction and experimental data of Table VII encouraging. Finally we compute total widths for D* and B* (for 
A H = 0.4 GeV): 

r(ZT°) = 38 KeV (83) 
T(D*+) = 62 KeV (84) 
T{B*° -> B°j) = 0.05 KeV (85) 



VI. CONCLUSIONS 



The increasing number of available data on heavy meson processes and the even more promising increase of data 
from forthcoming experiments demands theoretical predictions for these processes to be compared with experiment 
and to suggest subsequent lines of investigation. Calculating directly from the QCD lagrangian remains an extremely 
difficult task, in spite of the impressive success of recent lattice work. A most promising approach is the one based on 
heavy meson effective lagrangians, which incorporate the heavy quark symmetries and in addition the approximate 
chiral symmetry for light quarks. Although with increasing data such an approach will remain the best one beyond 
direct QCD calculations, at this stage it is made cumbersome by the large number of parameters that have to be fixed 
before obtaining predictions. In this note we have presented an intermediate approach, not as rigorous and general 
as that of the effective meson lagrangian, but which allows for a smaller number of input parameters. We start from 
an effective lagrangian at the level of mesons and of constituent quarks, and we then calculate the meson transition 
amplitudes by evaluating loops of heavy and light quarks. In this way we can compute the Isgur-Wise function, 
the form factors T\ji and T3/2, the leptonic decay constant F and F + , the coupling constants g and h relative to 
H — > Htt and H — > Sir processes respectively, the j3 coupling relative to the H — > processes. The agreement with 
data, when available, seems rather impressive. Additional data should be able to confirm some of the predictions of 
the model or suggest modifications. In conclusion the model presented here, based on meson and quark degrees of 
freedom, seems capable of incorporating a number of essential features of QCD and to provide a useful approach to 
calculate heavy meson transitions in terms of a limited number of physical parameters. 
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VII. APPENDIX 



We list here the integrals used in the paper. 



j AT /reg 



fc + A + ie 



/i 



16^/2 Jl/A: 

N c m 2 m 2 to 2 

iJVL f re9 d 4 k N r m 2 ^, m 2 m 2 



r eg ri 4 fc jV e m 2 m* 

16^ 4 7 (fc 2 - m 2 ) ~ 16^ 2 ( ' A 2 > 2 ' 



iN c f rcs , k 2 N c m 4 , 



16tt 4 J (fc 2 - to 2 ) 8tt 2 v ' A 2 ' [L 2 

iN c r cs d 4 k N c to 2 m 2 

2 16tt 4 y (fc 2 - to 2 ) 2 16^ 2 1 ' A 2 ' M 2 



J 3 (A) 



7 4 (A) = 



ijv„ r s d*k 



16tt 4 7 (fc 2 -TO 2 )(i;-fc + A + ie) 
AL f 1 ^ 2 



'l/A 2 

iiV c /- res d 4 fc 



16tt 4 7 (fc 2 -TO 2 ) 2 (wfc + A + ie) 



In these equations 



r(a,a; ,a; 1 ) = / di e _t 

J Xn 



is the generalized incomplete gamma function and erf is the error function. 
Let's introduce now: 

< a a ^ A! (l-x) + A 2 x 

a{x,A 1 ,A 2 ,w) - 



v/1 + 2 (w- 1) a; + 2 (1 - x 2 
Then the integrals used for the semilcptonic decays are as follows: 

iN c r s d 4 k 



-/ 

Jo 



16tt 4 7 (fc 2 - to 2 )(w • k + Ai + ie)(t/ • fc + A 2 + ie) 
1 



" l + 2a; 2 (l-cj) + 2a;(w-l) 
6 z^ 2 



167T3/ 2 ,/ 1/A 



/ ds a e - s (™ 2 -- 2 ) S -V2 (1 + erf^v^)) + 

./l/A 2 



6 



i/m 2 

l/A 2 



J 6 (Ai,A 2 ,w) 



16tt 2 

;al r os d A k 



iN c r cs 

16tt 4 y Jv~- 



16tt 4 7 (w • fc + Ai +ie)(t)' • fc + A 2 + ze) 



1\ I dx 



o 



l + 2a; 2 (l-w) + 2a;(w-l) 
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Nr. 



16tt 3 / 2 



dx 



l + 2x 2 (l-w) + 2x(w-l) 



^e-< m2 -^{a[l + erf(trv^)] ■ [1 + 2s(m 2 - a 2 
,/i/A2 L 



+ 2j-e" 

7T 



(95) 



S(Ai,A 2 ,w) - Ai/ 3 (A 2 ) +uj (h + A 2 / 3 (A 2 )) + Ai 2 7 5 (A 1 , A 2 ,w) 

T(Ai,A 2 ,w) = m 2 J B (Ai,A 2 ,w) + 7 6 (Ai,A 2 ,w) 

C/(Ai, A 2 ,w) = Ji + A 2 7 3 (A 2 ) + Ai J 3 (Ai) + A 2 Ax J B (Ai, A 2 , w) 



(96) 
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